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1.  Summary 


This  report  details  my  study  of  the  non-linear  plasma  oscillations  in  a  semiconductor  conduction 
channel  controlled  by  a  gate.  The  analysis  is  based  on  the  hydrodynamic  equations  derived  from 
the  Boltzmann  equation  for  a  parabolic  conduction  band.  The  hydrodynamic  approximation 
derived  here  includes  the  effects  of  viscosity,  finite  mobility,  and  temperature  gradients  in  the 
channel.  I  evaluated  the  electron-electron  collision  limited  mean  free  path  as  a  function  of 
temperature  and  gate  potential.  When  this  path  is  much  smaller  than  the  wavelength  of  the 
density  variations,  the  electron  gas  in  the  channel  can  be  treated  as  a  two-dimensional  fluid.  The 
flow  is  described  by  the  Navier-Stokes  equation  and  the  heat  conduction  equation.  The  quality  of 
the  plasma  resonance  is  limited  by  the  electron  mobility  and  the  viscosity  of  the  electron  fluid. 
Using  the  relaxation  time  approximation  I  evaluated  the  hydrodynamic  transport  coefficients, 
and  showed  that  at  low  temperatures  in  high  mobility  channels  the  quality  factor  of  the  plasma 
resonance  is  limited  by  the  viscosity.  The  hydrodynamic  nonlinearities  in  the  gated  channel 
subjected  to  a  time-harmonic  signal,  induce  a  constant  drain-to-source  voltage  with  a  resonant 
dependence  on  the  frequency.  This  effect  can  be  used  for  the  detection  of  electromagnetic 
radiation  by  the  Dyakonov-Shur  detector.  I  evaluated  the  detector  response  for  particular 
dimensions  and  values  of  mobility  at  which  plasma  resonances  are  in  the  terahertz  range.  I  also 
evaluated  the  quality  of  the  plasma  resonance  as  a  function  of  temperature,  and  showed  how  the 
viscosity  effects  limit  the  quality  factor  at  low  temperatures. 


2.  Introduction 


The  conduction  channel  of  a  semiconductor  Field  Effect  Transistor  (FET)  or  a  heterostructure 
High  Electron  Mobility  Transistor  (HEMT)  can  act  as  a  plasma  wave  resonator  for  density 
oscillations  in  quasi-two-dimensional  (2D)  electron  gas.  The  plasma  wave  here  refers  to  electron 
density  excitation,  possible  at  frequencies  significantly  higher  than  the  FET  cut-off  frequency  in 
a  short  channel  device.  The  hydrodynamic  model  predicts  a  resonance  response  to 
electromagnetic  radiation  at  the  plasma  oscillation  frequency  which  can  be  used  for  detection, 
mixing,  and  frequency  multiplication  in  the  terahertz  range  (1 ,  2).  In  particular,  the 
hydrodynamic  nonlinearities  produce  a  constant  source-to-drain  voltage  when  gate-to-channel 
voltage  has  a  time-harmonic  component.  In  the  Dyakonov-Shur  detector,  a  short  channel  HEMT 
is  used  for  the  resonant  tunable  detection  of  terahertz  radiation.  The  non-linear  plasma  response 
has  been  observed  in  InGaAs  (3,  4)  and  GaN  (5-3)  HEMTs,  in  the  frequency  range  from  0.2  to 
2.5  THz.  The  emission  of  radiation  is  also  possible  due  to  amplification  of  plasma  oscillations 
upon  reflection  of  plasma  waves  at  the  channel  lateral  boundaries,  which  has  been  shown  to  lead 
to  plasma  instability  under  certain  boundary  conditions  (9).  For  recent  reviews  of  the 


1 


experimental  efforts  in  detection  and  emission  of  THz  radiation  by  using  such  submicron 
channels  see  references  7  and  9.  The  results  show  that  short  channel  FETs  can  be  used  both  for 
tunable  resonant  and  broadband  detection  of  electromagnetic  radiation  in  THz  and  sub-THz 
range.  In  addition,  a  resonant  frequency  response  was  also  studied  in  a  few  GHz  range  using  a 
180  pm  long  channel  GaAs  HEMT  (10). 


The  plasma  waves  in  the  gated  two-dimensional  channels  have  linear  dispersion  law  (7) 

Mpi(q)  =  sq  =  q^ 


’(U.-U.) 


Ill 


(1) 


where  s  is  the  plasma  wave  velocity,  q  is  the  in-plane  wave  vector,  e  is  the  electron  charge,  Ug  is 
the  gate  voltage,  Uth  is  the  threshold  voltage  for  the  formation  of  the  conduction  channel,  and  m 
is  the  effective  electron  mass.  For  example,  if  gate  voltage  is  1  V  above  the  threshold,  the 

o 

plasma  wave  velocity  in  GaAs  and  GaN  channels  s  ~  10  cm/s,  usually  much  higher  than  the 
electron  drift  velocity.  Note  that  the  dispersion  law  in  the  gated  channel  is  different  from  the 
dispersion  law  in  the  un-gated  two-dimensional  electron  gas  (ungated  parts  of  FET)  with  sheet 
density  n0,  co(q)  =  (qe  n0/2m) 


In  a  FET  channel  with  length  L,  figure  1,  the  eigen-frequencies  of  the  plasma  standing  waves  are 
odd  multiples  of  the  fundamental  plasma  frequency,  given  by 


(2) 


Figure  1.  Schematic  structure  of  a  conduction  channel  in  a  FET.  The  gate  length  L  is  much  larger  than  the 
mean  free  path  of  electron-electron  collisions  A,ee  but  smaller  than  the  mean  free  path  of  electron- 
phonon  and  electron-impurity  scattering. 

If  the  mean  free  path  of  electron-electron  collisions  Xee  is  much  smaller  than  the  channel  length 
L,  the  plasma  transport  in  the  channel  can  be  studied  within  the  hydrodynamic  model  of  the 
electron  density  waves.  If  the  momentum  relaxation  time  tp,  determined  by  electron-phonon  and 
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electron-impurity  collisions  is  such  that  cdotp  »  1 ,  the  FET  can  operate  as  a  resonant  detector 
tunable  by  the  gate  voltage.  If  cootp<  1,  the  detector  response  will  be  broadband.  The  wavelength 
of  the  incoming  electromagnetic  radiation  in  the  terahertz  range  is  much  larger  than  the  device 
dimensions,  and  the  original  proposal  of  the  Dyakonov-Shur  detector  in  reference  1  included  an 
antenna  structure  in  order  to  collect  the  radiation.  In  a  design  with  a  slot  antenna,  an  ac  voltage 
Uac  is  induced  between  source  and  gate  contacts  as  shown  in  figure  2(a).  In  a  design  with  a  bow 
tie  antenna,  an  ac  current  Iac  is  induced  between  source  and  gate,  figure  2(b).  In  the  resent 
realizations  of  the  HEMT  based  detector,  references  3-7,  the  metal  pads  of  the  electrodes 
perform  the  role  of  the  antenna.  The  gate  voltage  Ug  is  tuned  to  control  the  electron  density  in  the 
gated  part  of  the  channel,  thus  allowing  the  resonant  plasma  frequency  to  vary  through  the 
frequency  of  the  incoming  signal. 


Figure  2.  Schematic  geometry  of  FET  operating  in  detector  mode  for  (a)  induced  ac  voltage  and  (b)  induced  ac 
current. 

The  analysis  based  on  hydrodynamic  model  of  electron  gas  in  the  conduction  channel  (7)  shows 
that  the  detector  response  AU,  the  constant  source-to-drain  voltage  induced  by  the  incoming  ac 
signal  of  amplitude  Ua  and  circular  frequency  to  in  the  configuration  shown  in  figure  2(a),  is 
given  by 
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(3) 


AU 

Un 


1  U  V 

^  a 

vU0y 


i(co). 


where  Uo=Ug-Uth.  In  high  mobility  short  channels  such  that  sxp/L  »1,  the  function  f(co)  has  a 
resonant  structure.  The  detector  response  in  the  above  equation  is  due  to  non-linear  plasma 
response  in  the  channel,  as  will  be  shown  in  some  detail  later  in  the  report.  The  propagation  of 
the  plasma  waves  in  the  channel  with  high  density  of  electrons  is  much  faster  than  electron  drift 
transport.  The  non-linear  properties  of  these  collective  excitations  allow  the  detector’s  response 
at  frequencies  considerably  higher  than  the  FET  cutoff  frequency,  determined  by  the  electrons 
ballistic  transport.  The  width  of  the  resonance  is  determined  by  the  momentum  relaxation  rate 
1/ip  that  limits  the  electron  mobility,  and  by  the  hydrodynamic  viscosity  in  the  channel.  In  the 
published  analysis  of  the  non-linear  plasma  waves  in  FET  channels  the  kinematic  viscosity 
coefficient  v  at  low  temperature  was  assumed  to  be  given  by  a  temperature  independent  value 
h/m,  about  15  cm  /s  for  GaAs  (1,  9). 

Below  I  will  derive  the  hydrodynamic  equations  from  the  quasi-classical  Boltzmann  equation 
using  relaxation  time  approximation  for  the  collision  integrals  and  obtain  the  temperature 
dependence  of  the  hydrodynamic  transport  coefficients,  i.e.,  viscosity  and  heat  conductivity. 
These  equations  will  be  derived  as  balance  equations  starting  with  a  non-equilibrium  Fermi- 
Dirac  distribution  as  a  zero  order  term  in  the  expansion  of  the  distribution  function  in  orders  of 
the  Knudsen  number,  the  ratio  of  the  mean  free  path  to  the  characteristic  length  of  density 
variations  (11).  The  first  order  corrections  to  the  pressure  tensor  will  be  obtained  using  the 
electron-electron  mean  free  path,  and  the  viscosity  tenns  will  be  obtained.  The  non-linear  terms 
in  the  hydrodynamic  equations  lead  to  non-linear  effects  in  the  plasma  wave  propagation  in  the 
channel.  In  particular,  a  constant  voltage  (zero  harmonic  tenn)  will  develop  in  the  plasma 
response  to  a  time-harmonic  perturbation.  I  will  apply  a  perturbation  analysis  to  the  evaluation  of 
this  non-linear  response,  evaluate  the  response  function  for  different  values  of  mobility  and 
viscosity,  and  estimate  the  detector  response.  It  will  be  shown  that  at  low  temperatures  in  high 
mobility  samples,  the  viscosity  is  the  main  physical  mechanism  that  limits  the  quality  of  the 
plasma  resonance.  I  will  include  the  heat  transfer  equation  in  the  model  and  evaluate  temperature 
variations  in  the  conduction  channel  for  relevant  values  of  physical  parameters. 


3.  Quasi-classical  Approximation 


Let  r  and  p  be  2D  position  and  momentum  and  v  is  2D  velocity  v  =  vFp/p.  The  quasi-classical 
Boltzmann  transport  equation  for  the  distribution  function  f(r,p,t)  of  electrons  in  the  2D 
conduction  channel  (z  =  0)  is 


— +  vVrf  +  eVrU-Vf 

a  r  r  p 


Vdty 


Vr 


(4) 
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where  U(r,z,t)  is  the  electric  potential  to  be  determined  self-consistently,  e  is  the  charge  of 
electron,  the  first  tenn  on  the  right  hand  side  is  the  collision  integral  due  to  electron-electron 
scattering,  the  second  term  is  the  collision  integral  due  to  electron-impurity  and  electron-phonon 
scatterings.  The  electron  density  n(r,t)  in  the  channel  is  obtained  from  the  distribution  function  as 

n(r,t)  =  gg  f  d2pf(r,p,t) ,  (5) 

(2ti h)  J 

2  is  spin  degeneracy  factor.  Electric  potential  is  found  from  the  Poisson  equation  (in 

V;U  +  ^  =  —  S(z)n(r),  (6) 

fe-  sss0 

—12 

where  ss  is  the  value  of  the  background  static  dielectric  constant  and  So  =  8.854x10  '  F/m.  If  the 
gate  contact  is  at  z  =  d,  the  boundary  condition  is  a  zero  tangential  field  at  z  =d.  In  the  limit  of 
the  infinitely  long  channel,  the  Poisson  equation  can  be  solved  by  Fourier  transformation, 

U(q,z)  =  J  d2re~'qrU(r,z)  (7) 


where  gs  = 
SI  units) 


with  the  result 


U(q,0)  = 


en(q) 

8s80q[l  +  coth(qd)] 


(8) 


In  the  long  wave-length  limit  qd  — >  0,  one  obtains  the  “gradual  channel”  approximation  (9)  U  = 
en/C,  where  capacitance  per  unit  area  C  =  ssSo/d. 


In  an  equilibrium  state  with  zero  drift  velocity,  f(r,p,t)  is  a  Fermi-Dirac  distribution  defined  by 
temperature  T  and  chemical  potential  We  assume  a  parabolic  conduction  band  and  obtain  from 
the  equation  (<9f/<9t)ee  =  0  (11): 


f°(p) 


1  +  exp 


^(p  -mu)2 
l  2kBT 


5 


(9) 


where  is  the  chemical  potential  C  divided  by  temperature,  =  (2/knT,  with  the  chemical  potential 
determined  by  the  normalization  equation  5  with  f  replaced  by  f° : 

q  =  kBTln(eEF/kBT-l).  (10) 

At  low  temperature  f°  —*■  O(p-pF)  where  Pf  is  Fermi  momentum.  We  also  define  Fermi  wave 
vector  kF  =  p i  / A  and  Fenni  energy  EF  =  c(T=0).  From  the  normalization  at  zero  temperature,  kF  = 

(27in)1/2  and 
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(11) 


„  nfi  n 

Ef= - . 

m 

It  is  important  to  note  the  applicability  of  the  quasi-classical  kinetics,  i.e.,  equation  4  at  low 
temperatutes.  The  spatial  condition  is  for  the  wave-length  of  the  applied  field  to  be  much  larger 
than  the  de  Broglie  wavelength  of  electrons:  1/k  »  ti!pv,  while  the  momentum  spread  should  be 

much  smaller  than  the  spread  of  the  distribution  f°(p):  Ti k  «  T/vF  where  vp=  pi  /m.  These 
conditions  and  the  time  domain  condition  are  equivalent  to  7/kvr  «  kBT,  tm «  EF .  Taking  k 
~  1/L,  the  spatial  condition  becomes 

kBT»2EF/LkF  (12) 


As  examples  consider  GaAs  and  GaN  channels  with  gate  to  channel  distance  d  =  35  mn.  For  the 
GaAs  channel  m/m0  =  0.067  and  ss  =  12.5.  The  gate  to  channel  capacitance  per  unit  area  C  = 
3.16x10  3  F/m2,  and  electron  density  n  =  UoC/e  where  Uo  =  Ug  -  Uth-  At  Uo  =  0.5  V  we  obtain  n 
~  lxlO12 1/cm2,  EF~  35  meV,  Vf  =  4.32  xlO7  cm/s.  If  L  =  0.2  pm,  equation  12  requires  T  »  16 
K.  For  the  GaN  channel  m/m0  =  0.2  and  ss  =8.9.  Then  C  =  2.25x10  3  F/m2.  At  Uo  =  1  V  we 
obtain  n  ~  1.4xl012  1/cm2,  Ep~  16.8  meV,  Vf  =  1.72xl07  cm/s.  If  F  =  0.2  pm,  equation  12 
requires  T  »  7  K.  Thus  the  validity  of  quasi-classical  theory  for  very  short  channels  requires 
temperatures  considerably  higher  than  10  K.  In  the  terahertz  range  f  ~  1  THz,  ho)  ~  4  meV,  and 

the  frequency  restriction  is  fulfilled. 


Considering  linear  response  to  a  hannonic  electrical  field  (77),  we  obtain  the  quasi-classical 
dielectric  function  in  a  two-dimensional  gated  channel: 


s(q,a>)  =  1  +  aK(q) 


00 

1  —  e  Ep/kBT  +  to  f  do 

'at;' 

0(co  -  qp/m)  0(qp/m  -  a>) 

0 

V®2  -  (qp/m)2  ^/(qp/m)2  -  ro2 

(13) 


where  0(co)  is  the  step  function, 


K(q)  = 


1 

qd[l  +  coth(qd)]  ’ 


e2md 

7is0ss/r 


(14) 

(15) 


At  temperatures  such  that  kBT  <  Ep/2,  cxp(— Ei  /k[jT )  is  small  and  the  integral  in  equation  13  can 
be  simplified  to  obtain  an  approximation 


s(q,(T>)  =  1  +  aK(q)  (l  -  e^^1 )  \  1 


a>0(a>-qvo)(l-GT)  /0(qvo -co)Gt 

I  2  2  2  l~2  2  2 

V®  -q  vo  Vo  vo  -® 


(16) 
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where 


GT(q,co) 


1-e  F  B  +e 


-Ef  /k„T 


(17) 


v0=vF 


1 


+  VIin(l_e-w) 
EF 


(18) 


The  plasmon  dispersion  is  obtained  from  the  equation  Re{s(k,copi)}  =  0.  Then,  if  the  factors  of 
exp(-EF/kBT)  are  neglected,  we  obtain 


®Pi(k)  = 


kvF _ 

1 

(l  +  l/aK(k))2 


(19) 


In  the  examples  considered  above,  a  is  much  larger  than  1,  namely  14.1  and  59.4  for  GaAs  and 
GaN  channels  respectively.  For  a  »1  and  kd  <  1,  equation  16  can  be  replaced  by 

»„i(k)»kvF^K(k).  (20) 

In  the  long  wave  length  limit,  kd  «1,  we  obtain  the  linear  dispersion  of  the  equation  1  copi(k)  = 
sk  with  s  =  (eUo/m)  ,  where  Uo  =  Ug-Uth.  It  follows  from  equation  19  that  copi  >  kvF  and 
Im  !c(k,c.0pi)}~  0.  Therefore,  in  the  quasi-classical  approximation  the  Landau  damping  (11),  a 
resonant  transfer  of  plasma  energy  to  single  particle  excitations,  is  small  in  the  gated  channel  for 
temperatures  below  EF/2. 


4.  Electron-electron  Scattering  Rate  and  Mean  Free  Path  in  Two- 
dimensional  Gated  Channels 


The  hydrodynamic  approximation  can  be  derived  from  the  kinetic  equation  assuming  a 
sufficiently  high  electron  density,  i.e.,  the  mean  free  path  of  the  electron-electron  (e-e)  collisions 
Xee  should  be  much  smaller  than  the  channel  length  L.  In  this  work  I  derive  the  Navier-Stokes 
equation  from  the  quasi-classical  Boltzmann  equation  (12,  13).  The  hydrodynamic  equations  can 
also  be  derived  by  statistical  averaging  of  the  Heisenberg  equations  of  motion  for  density  and 
momentum  operators  (14).  In  this  section  I  evaluate  the  mean  free  path  defined  by  electron- 
electron  collisions,  which  will  determine  the  limits  of  validity  of  the  hydrodynamic  equations. 
The  electron-electron  scattering  rate  will  be  evaluated  from  the  Fermi  Golden  Rule  as  done  for 
the  doped  channels  (15,  16),  using  quasi-classical  approximation  for  the  dielectric  function. 
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(21) 


1 


r«(P>o-) 


In 

ti 


k,q,cr' 


|y(q)|2 

e(q>(Ep-Ep+q)/a) 


x8(Ep.q  +  El.q-Ep-El) 


2  2 

where  index  a  refers  to  electron  spin,  E(q)  =  h"q“/2m,  v(q)  is  the  Coulomb  interaction  in  the 
gated  channel: 


v(q)  = 


sss0q  [l  +  coth(qd)]  ’ 


(22) 


and  the  exchange  terms  where  neglected  in  the  interaction  matrix  element.  The  inclusion  of  the 
exchange  terms  will  increase  the  scattering  rates  only  by  a  modest  amount  ( 1 7).  In  terms  of  the 
many-body  perturbation  theory,  the  first  term  in  the  square  brackets  in  equation  2 1  represents  the 
scattering  of  quasi-particles  and  the  second  term  represents  the  scattering  of  quasi-holes  (defined 
with  respect  to  the  Fermi  surface  of  the  electron  gas  in  the  conduction  band).  Summation  over  k 
and  a  in  equation  21  produces  the  imaginary  part  of  susceptibility  (15)  and  after  converting 
summation  to  integration,  we  obtain 


1 

'Ee(P) 


4n~h~ese0  0 


ldq7 

n  A 


1 


+  coth(qd) 


71 

c 

f F  -F  ^ 

d0 

coth 

p  p+q 

J 

0 

V  2ftkBT  J 

tanh 


( c-E  ^ 
2  p+q 

V  2f7kBT  J 


xlm - 

e(q,(Ep  -Ep+q)/ft) 


(23) 


where  C,  is  the  chemical  potential  given  in  equation  10. 

The  mean  free  path  for  an  electron  with  wave-vector  k  is  obtained  as  k(k)  =  VkXee(k)  where  Vk  = 
hk/m  and  xee  is  the  inverse  of  the  e-e  scattering  rate  in  equation  2 1 .  The  distribution  averaged 
value  of  the  mean  free  path,  after  summing  over  the  spin  directions,  is  defined  as 

—  1  1  00 

lTC*(PW  =  -  fdkkf“(k)A(k)  (24) 

*  zA.  Itn  o 

po 

where  p  denotes  the  momentum  while  k  denotes  the  wave-vector,  and  f° (k)  is  taken  at  zero  drift 
velocity.  The  mean  free  path  obtained  in  this  way  is  shown  in  figure  3  as  a  function  of 
temperature,  defined  in  units  of  the  Fermi  temperature  Tp  =  Ep/ks,  for  the  GaAs  channel.  The 
mean  free  path  for  the  GaN  channel  when  plotted  as  function  of  T/Tp,  is  similar.  It  follows  from 
figure  3  that  for  channel  of  length  L  =  0.2  pm,  the  hydrodynamic  approximation  for  transport  can 
be  used  for  T/Tp  >  0. 1 .  In  our  examples  the  Fermi  temperature  is  405  K  for  GaAs  and  195  K  for 
GaN  channels. 
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Temperature  in  units  of  the  Fermi  temperature 


Figure  3.  The  electron  mean  free  path  /,ee  determined  by 
the  electron-electron  scattering  is  shown  as  a 
function  of  temperature  for  the  GaAs  channel, 
when  the  gate  voltage  above  the  threshold  value 
is  U0=  0.5  V. 


5.  Derivation  of  the  Hydrodynamic  Equations  and  the  Transport  Coefficients 


The  hydrodynamic  equations  can  be  derived  from  the  Boltzmann  transport  equation  as  balance 
equations  for  conserved  quantities  (12).  The  hydrodynamic  variables  are  plasma  density  n(r,t), 
given  in  equation  5,  macroscopic  velocity  u(r,t),  and  internal  energy  S(r,t).  The  last  two  are 
obtained  by  an  ensemble  averaging  of  the  corresponding  microscopic  variables: 
u  =  (V),  8  =  ^m(v  -u)2/2^  .  For  a  parabolic  conduction  band  assumed  here,  the  microscopic 
velocity  and  momentum  are  linearly  related,  v  =  VpE(p)  =  p/m .  The  collision  integral  in 
equation  4  is  given  by  (11 ) 


Sf(  Pj) 


dt  j 


=  J</p2</p3</p4  {  ir(p3,p4,p1,p2)/(p3)/(p4)[l-/(p1)][l-/(p2)] 


(25) 


(Pi  ,  p2  ,  P3  ,  P4  )/(Pi  )/(P2  )[1  -  /(P3  )]  [1  -  /(P4  )]  }^(P3+P4-Pl-P2) 


2  n 


r(Pl,p2,p3,p4)  =  —  \M\  S(El+E2 -E, -E4)  , 


where  M  is  the  Coulomb  interaction  matrix  element.  The  average  of  a  microscopic  variable  Q  is 
defined  as 


Q  =  (Q) 


Jd2pf(r,p,t)0 

Jd2pf(r,p,t) 


(26) 
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The  balance  equations  for  the  density  (or  charge),  momentum,  and  energy  are  obtained  by 
setting  Q  =  1,  mv,  and  m(v-u)V2  respectively.  The  balance  equation  for  the  density  is  obtained 
in  the  fonn  of 


—  +  V-(nu)  =  0.  (27) 

dt 


The  electron-electron  collision  integral  does  not  contribute  because  of  charge,  momentum  and 
energy  conservation  in  equation  25  (12).  For  the  contribution  of  the  second  collision  integral  in 
equation  1, 1  will  use  the  relaxation  approximation  (14),  with  the  relaxation  time  x  being  an 
inverse  of  the  combined  rate  of  electron-impurity  and  electro-phonon  scatterings.  It  can  be 
obtained  from  measured  electron  mobility  in  the  channel.  This  gives  a  friction  tenn  in  the 
momentum  balance  equation: 

— +  (u-V)u  +  — VU  +  -  =  --V-fi  (28) 

dt  m  x  n 


where  the  pressure  tensor  was  defined: 

Hij  ^(^Vj-UjXVj-U;)). 


The  energy  balance  equation  takes  the  form 


n  —  +  nu  •  VS  +  V  •  q  +  TLA.  = 

dt  1  J 


v^tJc 


(29) 


(30) 


where  q  is  the  vector  heat  flux  defined  as 

q  =  ^<(v-„)(v-u)2), 


the  velocity  tensor  is  defined  as 


A, 


m  <3uj  fUj 
2  [fx,  +5xjy 


(31) 


(32) 


and  the  internal  energy  collision  integral  (<3£/St)c  is  related  to  the  energy  collision  integral 
(5W/0t)c  where  the  energy  W  =  S  +  mu2/2: 


c 

f  dW^ 

mu2 

Kdtj 

l  dt  J 

c  X 

(33) 


The  zero-order  approximation  to  the  balance  equations  is  obtained  by  replacing  the  distribution 
function  f  by  its  equilibrium  form  in  equation  9  where  the  density,  drift  velocity  and  temperature 
are  functions  of  position  and  time.  To  this  order  we  obtain  the  energy  in  terms  of  temperature 
and  density: 
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(34) 


£  = 


k.TF.Q;) 

F.(y 


where  the  Fermi  integrals  Fm  are  defined  as 


Fm(y)  =  Jdx 

0 


1  +  exp(x  -  y) 


(35) 


In  the  zero-order  approximation  the  heat  flux  vanishes,  q<0)  =  0,  while  the  pressure  tensor  is 
diagonal,  II(0)ij  =  6,jP,  where  the  scalar  pressure  is  given  by 


nS  nkBTF,(^) 
m  F0  (0 


(36) 


In  this  order  the  momentum  balance  equation  gives  the  hydrodynamic  Euler  equation  (18). 
Together  with  the  density  balance  equation,  it  was  used  to  predict  the  nonlinear  response  of  FET 
to  the  harmonic  perturbation  ( 1 ,  2,  9).  The  energy  balance  equation  can  be  put  in  form  of 
equation  for  electron  temperature.  Then,  use  the  energy  units  and  define  9(r,t)  =  kBT(r,t).  Define 
the  heat  capacity  as  cv  =  (<3£/59)n  =  (<9S/3T)n/kB  and  evaluate  it  as  the  heat  capacity  of  the  two- 
dimensional  ideal  Fenni  gas: 


29F,(g)  Ef 
Ef  e(l-e-ve) 


(37) 


where  EF  is  a  function  of  n,  see  equation  1 1 .  Now,  writing 


V£  = 


(d£) 

Vn  + 

f  58^1 

v5ny 

0 

laej 

ve 


and  using  equations  27  and  34  I  obtain  the  equation  for  the  electron  temperature  distribution. 
Thus,  to  zero-order  we  obtain  the  Euler  equation  and  the  temperature  distribution  equation: 

— +  (u-V)u  +  — VU  +  -VP  +  -  =  0 
dt  m  n  x 


50  _  ,  if  5W3  mu2 

—  +  V-(0u  = —  - - 

dt  cl  dt  L  c„x 


(38) 

(39) 


From  equations  10  and  36  the  hydrostatic  pressure  can  be  obtained  as 

P  jim 

ntr 


(40) 


The  effects  of  viscosity  and  thermal  conductivity  appear  in  the  next  order  of  expansion  in 
Knudsen  number  (12).  Correspondingly,  write  the  distribution  function  as 
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f(r  ,v  ,t)  =  f 0  (r  ,v,t)  +  f  < 1 '  (r ,  v,t) 


(41) 


A  systematic  expansion  of  the  classical  distribution  function  for  a  parabolic  band  exists  in  the 
form  of  the  Chapman-Enskog  expansion  (12).  Here  I  will  use  a  relaxation  time  approximation 
for  the  first  order  correction: 


'  df  \  ~  f(1) 

UtJee  Ue 


(42) 


Then,  from  the  Boltzmann  equation  with  only  the  electron-electron  scattering  included,  one 
obtains  to  this  order 


fd)  _ 


— 


— +  vVrf°+eVrU-V  f° 
St  p 


\ 

) 


(43) 


with  f°  =  f  °(n(r,t),u(r,t),0(r,t))  ,  where  n(r,t)  and  u(r,t)  are  solutions  of  the  zero  order 
hydrodynamic  equations.  Without  the  relaxation  time  approximation,  one  has  to  solve  the 
linearized  Boltzmann  equation  in  order  to  obtain  the  transport  coefficients  in  the  first  order  of  the 
expansion  in  Knudsen  number  (11).  In  case  of  fennions  in  ungated  two-dimensional  channel,  an 
analytical  solution  is  possible  in  the  low  temperature  limit  (19).  A  much  simpler  approximation 
in  equation  42  allows  one  to  obtain  an  estimate  of  the  hydrodynamic  transport  coefficients  in  the 
wider  temperature  range  while  giving  the  correct  temperature  dependence.  The  derivation  is 
outlined  below.  After  taking  the  partial  derivatives,  the  first  order  contributions  are  evaluated  at 
u  =  0  and  we  obtain 


f"’=-T„(p)[f0(l-f") 


m  2 

— -w‘  - 
202 


2F, 


'F  j 


1  ( 

W-V0  +  — 

0 


W,W  ■ - 8::W 

1  1  2  1 


(44) 


where  w  =  v  -  u.  Using  this  in  the  evaluation  of  the  ensemble  averages  in  equations  29  and  31, 
one  obtains  the  pressure  tensor  and  the  heat  flow  to  the  first  order.  The  pressure  tensor  can  be 
written  as  a  sum  of  a  unit  tensor  and  a  traceless  symmetric  tensor  (18),  n  i,  =  5ijP  +  n(1)ij.  Let  us 
define  the  following  integrals  over  the  momentum  dependent  relaxation  time: 

T?  00 

>  m(T)  =  — L  Jdkk-x„(k)f"(k)[l  -  f°(k)]  (45) 


and  obtain  the  first  order  term  in  the  pressure  tensor  in  the  fonn 


nij1)  =  ■ 


2vn 


m 


A 

Aij~y5uv'u 


V  ^  J 

where  the  kinematic  viscosity  coefficient  v  is  given  by 


(46) 
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(47) 


AIl 

m  T 


J5(T). 


Notice  that  the  second  viscosity  (18)  does  not  appear  in  the  pressure  tensor  in  this  order,  as  also 
is  the  case  in  the  classical  statistics  with  a  parabolic  energy  band  (11). 


The  heat  flow  is  obtained  from  equations  3 1  and  44,  and  to  the  first  order  it  is  proportional  to  the 
temperature  gradient,  as  also  is  the  case  in  the  classical  fluid: 

q  =  -  kV0  ,  (48) 


where  the  coefficient  of  thermal  conductivity  k  is  given  by 


K 

n 


2 n 

m 


J  T 


(49) 


Using  equations  46  and  48,  one  obtains  the  equations  of  the  viscous  hydrodynamic  model  for  the 
electron  fluid  in  the  gated  channel  as  the  two-dimensional  density  balance  equation,  the  Navier- 
Stokes  equation,  and  the  heat  equation.  Omitting  the  gradients  of  the  transport  coefficients,  the 
equations  are: 

~  +  V '  (nu)  =  0  (50) 

dt 


— +  (u-V)u  +  — VU  +  -VP  +  --vV2u  =  0 
dt  m  n  x 


— +  V-(0u) 
dt  {  ’  c 


2Lv20. 


mv 


2c 


3m  dm 
— L  +  — - 

OX  OX; 


■  8;;  — — -k 
3x, 


f  8W  ^ 

l  5t  J 

mu 


c,  x 


(51) 

(52) 


where  we  defined  the  coefficient  y  =  K/n,  which  has  the  same  units  as  kinematic  viscosity.  This 
coefficient  divided  by  the  heat  capacity,  y/cv,  is  the  thermo  metric  conductivity  (18).  The 
coefficients  of  viscosity  and  thennal  conductivity  for  GaAs  and  GaN  channels  are  shown  in 
figure  4  as  functions  of  temperature,  T/TF. 
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Figure  4.  Coefficients  of  kinematic  viscosity  (left  panel)  and  thermal  conductivity  divided  by  density  (right 
panel)  are  shown  as  functions  of  temperature  for  the  (a)  GaAs  and  (b)  GaN  channels. 


In  the  low  temperature  limit  the  electron-electron  collision  time  Tee(p)  in  equation  45  can  be 
replaced  by  its  value  at  the  Fermi  surface,  xee(PF)  {15,  16)  and  an  explicit  temperature  dependence 
can  be  obtained  for  the  transport  coefficients: 


2h  T2 

v(T  «  TF)  ~ 


ran  T2  ln(2TF/T) 


k(T«Tf)  4nti  Tf  1 


n 


3m  T  ln(2TF/T) 


(53) 

(54) 


Notice  that  the  temperature  dependence  is  same  as  one  would  obtain  without  the  relaxation  time 
approximation  by  solution  of  the  linearized  Boltzmann  equation  {19). 

It  is  clear  from  figure  4  and  equation  53  that  at  moderately  low  temperatures  the  viscosity  can  be 
much  larger  than  the  value  assumed  in  the  published  studies  of  FET  based  plasma  resonator,  a 
constant  equal  h/m  {1-9).  As  shown  below,  at  the  temperatures  at  which  many  of  the 
experimental  studies  of  Dyakonov-Shur  detector  were  performed  {3-7,  10)  the  viscosity  can  be  a 
major  limiting  mechanism  for  the  quality  of  the  plasma  resonance. 


In  order  to  compare  the  relative  contributions  of  the  self-consistent  field  and  the  pressure 
gradient  terms  in  the  Navier-Stokes  equation  51,  assume  a  uniform  flow  and  use  the  gradual 
channel  approximation,  U=  en/C,  and  then  obtain: 


(l/n)VP  ^  ntrC  ]  |  kBT  F,  F, 
(e/m)VU  ~  me2  l^F^  (l-e^)F0 


(55) 


This  ratio  generally  will  be  small  if  T  <  Tp.  Consider  for  example  a  GaAs  channel,  d  =  35  nm,  at 
Ug  -  Uth  =  1  V.  Then  the  low  temperature  limit  of  the  right  hand  side  in  equation  55  is  0.075.  The 
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linear  approximation  to  the  hydrodynamic  equations  has  solutions  in  the  form  of  plasma  waves 

/ 2  2 

propagating  with  the  speed spl  =  ^so  +  s«  ,  where  So  is  obtained  from  the  electric  field  term  in 

the  momentum  balance  equation  and  is  given  in  equation  20,  and  sa  is  the  speed  of  sound 
obtained  from  the  pressure  gradient  term.  At  the  electron  density  in  the  channel  with  gate 
voltages  of  0.5  =  10  V,  sa  «  so" ,  one  can  neglect  the  pressure  gradient  term  compared  to  the 
electric  field  tenn  (20). 

The  relaxation  time  x  in  the  friction  term  in  equation  50  can  be  related  to  the  mobility  q  in  the 
channel,  q  =  ex/m.  For  the  collision  tenn  in  equation  52,  (fW/ft)c,  I  will  use  a  relaxation 
approximation  (13),  (0  -  9o)/xe  where  9o  is  the  lattice  temperature,  and  the  energy  relaxation  time 
Tg  can  be  calculated  from  the  electron-phonon  scattering  rate  (14). 


6.  Nonlinear  Response  of  the  Channel  Confined  Electron  Plasma  to  a  Time- 
harmonic  Signal 


As  a  simple  application  of  the  hydrodynamic  equations  in  this  section,  I  will  calculate  the 
response  function  of  the  Dyakonov-Shur  detector  based  on  the  FET  conduction  channel  using  the 
equivalent  circuit  in  figure  2(a).  In  doing  this,  I  will  generalize  the  treatment  in  reference  1  to 
include  the  effects  of  viscosity,  and  also  evaluate  the  electron  temperature  distribution  in  the 
channel.  The  boundary  conditions  are  U(0,t)  =  U0  +  Uacos(oot)  at  the  source  side  of  the  channel 
and  zero  current  at  the  drain  side,  j(L,t)  =  0.  In  the  detector  mode,  the  time-harmonic  part  is 
induced  by  the  external  radiation  (1-7).  Here  the  analysis  is  restricted  to  a  one-dimensional  flow 
that  is  unifonn  in  the  transverse  channel  direction: 


fn  f  (nu) 

—  +  — — -  =  0 
ft  dx 


(56) 


fu  <3u  e  fU  1  fP  u  f 2 u  . 

—  +  u  —  + - + - +  —  v — 7  =  0 

ft  5x  mfx  n  fx  x  fx“ 


f0  f(0u)  x  S20  mv 

6fuA 

(2_J_ 

6  fW> 

mu2 

ft  fx  cv  fx2  cv 

v  fx  J 

Cv 

l  St  J 

c  CA 

(57) 

(58) 


In  the  narrow  channel,  approximation  n  =  CU/e.  Then,  the  boundary  condition  at  the  source 
contact  becomes  n(0,t)  =  n0  +  nacos(cot).  As  agued  above,  the  pressure  gradient  term  in  equation 
57  can  be  omitted  as  compared  to  the  electric  field  term.  Then,  if  one  assumes  that  the 
temperature  variation  is  small  enough  to  ignore  the  space-time  variation  of  the  viscosity,  the 
Navier-Stokes  equation  effectively  decouples  from  the  heat  equation,  equation  58. 

For  convenience  use  dimensionless  variables  and  parameters,  defined  as 
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(59) 


t  =  ts0/L,  x  =  x/L,  n  =  n/n0,u  =  u/s0,x  =  xs0/L,v  =  v/(s0L),  %  =  x/Ls0,  0  =  0/mso 


where  So  is  the  equilibrium  value  of  the  plasma  wave  velocity.  The  hydrodynamic  equations 
derived  in  the  previous  section  can  be  written  in  the  dimensionless  form, 

dh  d(nv)  . 
dt  dx 

i  .i  ffi  u  .  52u  . 

— +  u  —  +  —  +  —  v — -  =  0, 
dt  dx  dx  x  dx2 


50  5(0  u)  x  s2e 

dt  dx  c„  dx2 


-  v 


5u 

\dxj 


To  C„X 


(60) 

(61) 

(62) 


where  we  defined  xe  =  cvxe  .  When  the  solution  u(x,t)  of  the  equations  60  and  61  is  found,  it  can 
be  used  in  equation  62  to  obtain  an  inhomogeneous  heat  equation  which  then  can  be  solved  to 
obtain  the  electron  temperature  distribution  0(x,t)  in  the  channel. 

The  boundary  conditions  are 


n(0,t)  =  l  +  nacos(eot),  u(l,t)  =  0. 

(63) 

Expand  the  solutions  of  equations  59  and  60  in  time  harmonics: 

n-l  =  nm  +  nm+nl2)+... 

(64) 

d  =  u(0)  +  u(1)+u(2)+... 

(65) 

where  n(m),v(m)  oc  exp(/mcot ),  m=0,l,2. ...  The  hydrodynamic  nonlinearities  result  in  coupling 
of  different  harmonics.  Assuming  that  the  input  signal  is  small,  na  « 1  and  u  « 1 ,  apply  the 
small  signal  analysis.  This  obtains  the  following  hierarchy  of  orders  of  expansion: 

n(1,,u(1)  =  0(na);  n<V0)  =  0(fi2); u(0\u(2)  =  0(fi2) .  (66) 


Separating  the  orders  and  integrating  over  the  time  in  the  second  order  equations,  obtain 

5nU)  5u(1) 


5t 

.0) 


5u  5n'"  u 

— z — i - 1 - v - - 

dt  dx  x  dx~ 


dx 

x 


=  0 

_  52u(1) 


=  0 


_5 

5x 


^(u(0)  +(n(1)u(1)))  =  0 


(67) 

(68) 

(69) 
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where  the  angular  brackets  denote  time  average  over  the  period  2tc/co.  Boundary  conditions  for 
these  equations  are 

n(1|(0,t)  =  nacos(cot),  n,0)(x  =  0)  =  ul0)(x  =  1)  =  u(1)(l,t)  =  0. 

The  induced  constant  source  to  drain  voltage  will  be  found  as 

A U  =  (e/C))n0An(0)  =(e/C))n0  [n(0)(x  =  1  )-n(0\x  =  0)]  (71) 

From  equations  69  and  70,  obtain 

An(0)  =-(u(1)2(0))  +  ^Jdx(n(1V1,)-v(n(1)(0)u(1,(0))  (72) 

2  t  „ 


where  the  arguments  refer  to  x.  The  first  order  functions  n(1)  and  u(1)  are  solutions  of  equations  67 
and  68  which  give  linear  wave  equation,  with  solutions  proportional  to  exp(z'kx-uot).  The 
dispersion  relation  is 


of  +  mix  -  k2  (s02  -  moj  =  0 . 

(73) 

The  solutions  of  equations  67  and  68  are 

n(1)  =Re|(cie'k°x  +C2e’/k°x)e"!<yt} 

(74) 

u(1)  =  Re  j(C1e*k°x  -C2e‘&°x )  e"'®‘ } 

(75) 

where  k0(co)=ki  +  zk2  is  solution  of  equation  73  and 


Here  I  defined  ii  =  v/s02.  Constants  Ci  and  C2  are  found  from  the  boundary  conditions: 


r 


l  +  exp(±2zk0L) 


Then  use  equations  74  and  75  in  equation  72,  and  use  equation  71  to  find  constant  source  to 
drain  voltage  AU  induced  by  the  harmonic  input  signal.  This  can  be  written  as  in  equation  3: 
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AU 

Un 


^  u  v 

^  a 

vU0y 


where  the  detector  response  function  f(oo)  is  given  in  terms  of  ko(oo).  Let  us  define  the  following 
combinations  of  hyperbolic  and  trigonometric  functions  of  kjL  and  k2L: 

Sl  -  cosh(2£2L)  -  cos(2^L) 


S2  =  —  [cosh(2k2Z)  - 1]  +  —  [1  -  cos(2k,Z)] 
k  2  k  ■ 


53  =  cosh(2 k2L)  +  k2  cos(2k\L) 

54  -  cosh(2k2L)  +  cos(2klL) . 


In  terms  of  these  the  detector  response  function  is  given  by 


fM 


(ox)2  S[  coxS2  2VWXS, 

1 

l 

cx 

<N 

ex 

_ 1 

k0L 

(78) 


where  v  and  x  are  defined  in  equation  59.  Depending  on  the  value  of  the  friction  parameter  i,  the 
detector  response  will  be  resonant  or  broadband. 

Next,  consider  the  behavior  of  the  response  function  in  the  case  of  the  GaAs  channel  of  length  L 
=  0.5  pm,  gate  insulator  thickness  d  =  35  mn,  and  set  Uo  =  0.5  V  as  in  the  first  example  in  section 
III.  The  temperature  dependence  of  the  friction  parameter  at  low  to  moderate  temperatures  is 
determined  by  the  electron-acoustic  phonon  scattering  rate,  and  the  total  rate  is  given  by 

-  =  —  +  aT  (79) 

x  xn 


where  the  zero  temperature  value  1/xo  is  due  to  scattering  from  impurities  and  defects.  We  set 

5  2 

here  xo  =  10  ps  (the  corresponding  low  temperature  mobility  is  about  2.7x10'  cm  /Vs).  For  GaAs 

Q  #  #  0 

heterostructure  we  obtain  (21,  22)  a  =  7. 1 x  10  1/Ks.  This  detennines  the  friction  parameter  i  as  a 
function  of  temperature,  while  the  temperature  dependence  of  the  viscosity  coefficient  is  found 
from  equation  47  and  figure  4.  The  resonant  frequency /o  =  t0o/2n=  so/4L  =  0.57  THz.  In  figure  5 
the  resonant  response  is  shown  at  two  values  of  T/TF,  0.1  and  0.15,  corresponding  to  T  =  40.6  K 
and  60.9  K.  The  corresponding  values  of  the  friction  parameter  Sox/L  are  17.7  and  15.9  and  of  the 
kinematic  viscosity  coefficient  260  and  135  cm  /s.  The  higher  resonances  are  strongly 
diminished  by  the  viscosity  effect. 
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Figure  5.  Response  function  of  the  GaAs  channel  based  detector  at 
two  different  temperatures:  T/TF  =  0.1  (solid  line)  and 
T/Tf  =  0.15  (broken  line). 

When  the  friction  parameter  xso/L  decreases  below  0.5,  the  response  becomes  broadband, 
corresponding  to  overdamped  plasma  oscillations  in  the  conduction  channel,  as  shown  in  figure  6 
for  tso/L  =  0.3  and  0.2  with  v/soL  =  0.001.  For  a  long  channel,  as  xso/L  decreases  the  plasma 
waves  decay  before  ever  reaching  the  drain  contact. 

3.0 
2.5 
2.0 
3^1-5 
1.0 
0.5 
0.0 

COT 


Figure  6.  Response  function  in  the  case  of  strongly  damped  plasma 
oscillations,  giving  a  broadband  detector  response,  for  two 
different  values  of  friction  parameter,  tso/L  =  0.3  and  0.2  with 
v/s0L  =  0.001. 
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Next,  turn  to  the  evaluation  of  the  electron  temperature  distribution  in  the  channel  to  0(na).  To 

2  7 

that  order  the  u~  and  (Su/fixfi  terms  in  equation  62  can  be  omitted.  Define  A0  =  0  -  0o  and 

®(x,t)  =  -0ofi  (80) 

ox 


where  0o  refers  to  lattice  temperature  and  u(x,t)  is  given  in  equation  74.  Then  obtain  an  equation 
for  A0: 


fiAO  _„,fi2A0 
fit  ^  fix2 


+  ®(x,t) 


(81) 


where  %'  =  yj cv .  For  the  boundary  conditions,  assume  lattice  temperature  at  the  source  and  drain 
contacts  and  assume  that  the  electron  temperature  was  initially  set  to  the  lattice  value: 

A0(O,t)  =  A0(  1  ,t)  =  0,  A0(x,0)  =  0 .  (82) 

Substitution  A0(x,t)  =  f(x,t)exp(- t/xe)  transforms  equation  80  into 


fif  fi2f 

=x'^rr  +  <l)(x,t)exp(-t/xe),  (83) 


which  is  the  inhomogeneous  heat  equation  with  Dirichlet  boundary  conditions  (23).  The  Green’s 
function  for  this  equation  can  be  evaluated  by  Fourier  expansion  in  the  interval  0<  x  <1,  and 
restoring  the  dimensional  units  we  obtain 


AT 


m=l 


[ exp  (-/cot)  —  exp(  -  a mt) ] 


k0L7imsin(7imx/L)  j^l  +  (-l)m  1sec(k0L) 


(o)L/s0  +/am)(7i2m2  -k2L2) 


(84) 


where  am 


_E2x' 

Ls 


L 

s<Te 


As  an  example,  consider  a  signal  with  the  amplitude  equal  to  0.1  of  the  gate  voltage  applied  to  a 
GaAs  channel  with  sox/L  =  2  (corresponding  to  mobility  of  54000  cnr/Vs),  and  we  set  SoXe/L  =  2, 
v/soL  =  0.04,  xAcySoL)  =  0.16.  The  deviation  of  the  electron  temperature  T(x,t)  from  the  lattice 
temperature  To  is  shown  in  figure  7  as  a  function  of  the  position  along  the  channel  at  times  tso/L 
=  1  and  2,  evaluated  at  the  resonance  frequency. 
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Figure  7.  An  example  of  electron  temperature  distribution  in  the 
conduction  channel  with  s0t/L  =  2,  at  two  different  times 
(a)  tSo/L  =  1  and  (b)  tS(/L  =  2.  Both  distributions  are 
evaluated  at  the  resonance  frequency  co0,  with  the  amplitude 
of  the  applied  signal  equal  to  0. 1U0. 


As  another  example,  consider  the  case  corresponding  to  the  solid  line  in  figure  5,  i.e.,  GaAs 
channel  at  lattice  temperature  T  =  0.1  Tp,  when  the  signal  of  amplitude  0.0  lUo  and  to  =  too  is 
applied.  The  resulting  temperature  distribution  is  shown  in  figure  8,  at  two  times,  (a)  tso/L  =  1 
and  (b)  tso/L  =  2. 


x/L 


Figure  8.  Electron  temperature  distribution  in  the  GaAs  conduction 

channel  at  T0  =  0.1TF  at  two  different  times,  (a)  tso/L  =  1  and 
(b)  tso/L  =  2  in.  Both  distributions  are  evaluated  at  the  resonance 
frequency  co0,  with  the  amplitude  of  the  applied  signal  equal  to  0.0 1U0. 
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Note  that  in  these  examples  the  variation  of  the  electron  temperature  in  the  channel  is  relatively 
small  when  the  applied  signal  is  small.  This  justifies  the  approximation  of  uniform  electron 
temperature  in  the  solution  of  the  Navier-Stokes  equation.  At  much  higher  mobility  and  when  the 
applied  signal  is  not  small  as  compared  to  the  gate  voltage,  the  temperature  dependence  of  the 
transport  coefficients  may  have  to  be  included  explicitly  in  equations  60-62  which  then  would 
have  to  be  solved  numerically. 

One  can  estimate  the  sensitivity  of  the  detector  following  reference  1 .  If  the  electro-magnetic 
radiation  is  coupled  to  device  by  means  of  dipole  antenna,  the  detector  response  is  given  by  R  = 
AU/(1S),  where  1  is  the  radiation  energy  density  in  the  incoming  signal,  AU  is  given  in  equation 
3,  and  the  antenna  aperture  S  =  a”G/4tc.  Here  A  is  the  wavelength  and  G  ~  1.5  is  the  antenna  gain. 
We  then  obtain  an  estimate  R  ~  ( 1 00/XJ0)f(co)  V/W.  For  the  example  shown  in  figure  5,  R  ~ 

104  V/W. 


The  predicted  response  is  smaller  than  was  originally  predicted  for  semiconductor  HEMT  device 
at  low  temperature  (7),  due  to  significant  increase  of  viscosity  coefficient  with  the  decreasing 
temperature.  Most  experiments  with  the  Dyakonov-Shur  plasma  resonators  have  been  done  at 
low  temperature  in  order  to  reduce  electron-phonon  scattering  rate  1/x  and  achieve 
correspondingly  high  value  of  mobility,  see  references  1-10.  It  follows  from  the  present 
calculation  that  at  very  low  temperatures  the  entropy  production  associated  with  viscosity 
negates  the  advantage  of  increasing  carrier  mobility  with  lower  temperature.  One  can  see  it  in 
considering  the  quality  of  the  plasma  resonance  at  low  temperatures.  Near  the  resonant  frequency 
®o  the  function  f(co)  is  approximately  a  Lorentzian  with  the  width  Am  detennined  by  the  friction 
and  viscosity.  Define  the  quality  factor  Q  of  the  resonance  at  m  =  Mo  as  Q  =  Mq/Am.  Then 


1  2L  27iv 

-  — - 1 - . 

Q  7ts0i  s0L 


(85) 


The  temperature  dependence  of  v  is  given  in  equation  47  and  figure  4,  and  the  temperature 
dependence  of  x  at  low  temperature  is  given  in  equation  79.  The  linear  coefficient  a  in  equation 
79  can  be  deduced  from  the  temperature  dependence  of  electron  mobility  in  GaAs  {21,  22)  and 
GaN  {24)  channels,  and  we  obtain  values  of  7. 1  x  10'  1/Ks  and  22x10'  1/Ks,  respectively.  Then 
the  equation  85  can  be  rewritten  as 

1  2L  2LaT  2nv{T) 

—  = - + - + - — .  (8 

Q  71SqXq  71Sq  s0L 


With  the  same  parameters  as  in  figure  4  and  Xo  =  10  ps,  the  resonator  quality  factor  Q  is  shown  in 
figure  9  as  a  function  of  temperature  at  low  temperatures  for  GaAs  and  GaN  submicron  channels, 
L  =  0.5pm,  with  the  plasma  resonance  at  approximately  0.5  THz. 
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Temperatute  in  units  of  the  Fermi  temperature 


Figure  9.  Plasma  resonator  quality  factor  is  shown  as 

function  of  temperature  for  (a)  GaAs  and  (b)  GaN 
channels  with  length  L  =  0.5  pm.  The  resonance  is 
at  f0  ~  0.5  THz. 


One  can  see  that  the  quality  factor  is  limited  to  values  below  Q  =  10,  which  agrees  with  the 
experimental  studies  ( 4 )  in  which  the  loss  to  transversely  propagating  waves  (oblique  plasmons 
(25,  26))  was  eliminated  by  having  a  multi-channel  configuration. 


Next,  consider  resent  low  temperature  experiments  (10)  with  the  long  channel  GaAs  detector, 
with  L  =  180  pm,  d  =  190  nm,  and  the  measured  electron  density  no  =  1.3x10  1/cm  . 
Correspondingly,  the  plasma  resonance  is  at  fo  ~  1 .2  GHz.  By  fitting  the  measured  values  of 
mobility  as  1/p  =  1/po  +  aT  and  relating  x  to  p  by  x  =  mp/e,  we  find  that  the  scattering  rate  is 

o 

obtained  from  equation  79  with  xo  =  0.539  ns  and  a  =  5.45x10  1/Ks.  The  viscosity  coefficient 
can  be  calculated  from  equation  53  as 


v(T)  «  1 1  (cm2/s): 


f  t  3 
_!jl 

v  T  y 


ln(2TF/T) 


where  Tp  =  54  K.  The  quality  factor  evaluated  from  equation  85  is  shown  in  figure  10  and  the 
results  are  in  reasonable  agreement  with  the  resonance  widths  in  reference  10. 
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Figure  10.  Plasma  resonator  quality  factor  is  shown  as  function  of 
temperature  for  a  GaAs  detector  with  channel  length 
L  =  180  pm.  The  resonance  is  at  f0  ~1.2  GFlz. 

Thus,  as  in  the  case  of  short  channel  plasma  resonators  discussed  above,  here  too  the  viscosity 
effects  limit  the  quality  of  the  resonator  at  low  temperatures,  reducing  the  advantage  of  higher 
mobility  achieved  by  lowering  the  temperature. 


7.  Conclusions 


In  this  work  I  considered  the  non-linear  plasma  transport  in  gated  conduction  channels.  The 
hydrodynamic  Navier-Stokes  equation  was  derived  from  the  quasi-classical  Boltzmann  equation, 
and  the  applicability  of  hydrodynamic  model  was  evaluated.  The  non-linear  convection  tenn  in 
the  Navier-Stokes  equation  and  the  non-linearity  in  the  current  tenn  in  the  mass  balance  equation 
contribute  to  non-linear  effects  in  the  wave  propagation  in  the  electron  fluid.  I  evaluated  the 
temperature  dependent  hydrodynamic  transport  coefficients  in  gated  two-dimensional  conduction 
channels.  I  applied  these  equations  to  show  that  a  constant  voltage  develops  in  the  plasma 
response  to  a  time-hannonic  perturbation,  and  evaluated  the  response  function  for  different 
values  of  mobility  and  viscosity  and  we  then  estimated  the  detector  response.  The  treatment  here 
extended  the  theory  of  Dyakonov-Shur  plasma  resonator  and  detector  (7)  to  account  for  the 
temperature  dependence  of  viscosity,  and  also  included  the  energy  balance  equation  into  the 
analysis.  Depending  on  the  length  of  the  gated  channel,  the  detector  response  can  be  in  the  range 
from  a  few  GHz  to  a  few  THz,  and  is  further  tunable  by  the  gate  voltage.  Depending  on  electron 
mobility  in  the  channel  the  response  can  be  resonant  or  broadband.  I  showed  that  the  quality  of 
the  plasma  resonance  in  the  high  mobility  channels  at  low  temperature  is  limited  by  the  effects  of 
viscosity.  Though  the  evaluation  of  the  response  was  confined  to  the  one-dimensional  uniform 
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flow,  the  hydrodynamic  equations  derived  here  can  also  be  applied  to  the  two-dimensional  wave 
propagation,  and  the  effects  of  the  oblique  modes  ( 26)  can  also  be  included. 
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No.  of 

Copies  Organization 

1  ADMNSTR 

ELECT  DEFNS  TECHL  INFO  CTR 
ATTN  DTICOCP 

8725  JOHN  J  KINGMAN  RD  STE  0944 
FT  BELVOIR  VA  22060-6218 

1  CD  OFC  OF  THE  SECY  OF  DEFNS 
ATTN  ODDRE  (R&AT) 

THE  PENTAGON 
WASHINGTON  DC  20301-3080 

1  US  ARMY  RSRCH  DEV  AND  ENGRG 
CMND 

ARMAMENT  RSRCH  DEV  &  ENGRG 
CTR 

ARMAMENT  ENGRG  &  TECHNLGY 
CTR 

ATTN  AMSRD  AAR  AEF  T  J  MATTS 
BLDG  305 

ABERDEEN  PROVING  GROUND  MD 
21005-5001 

1  PM  TIMS,  PROFILER  (MMS-P) 
AN/TMQ-52 
ATTN  B  GRIFFIES 
BUILDING  563 
FT  MONMOUTH  NJ  07703 

1  US  ARMY  INFO  SYS  ENGRG  CMND 
ATTN  AMSEL  IE  TD  A  RIVERA 
FT  HUACHUCA  AZ  85613-5300 

1  COMMANDER 

US  ARMY  RDECOM 
ATTN  AMSRD  AMR 
WC  MCCORKLE 
5400  FOWLER  RD 

REDSTONE  ARSENAL  AL  35898-5000 

1  US  GOVERNMENT  PRINT  OFF 

DEPOSITORY  RECEIVING  SECTION 
ATTN  MAIL  STOP  IDAD  J  TATE 
732  NORTH  CAPITOL  ST  NW 
WASHINGTON  DC  20402 

1  US  ARMY  RSRCH  LAB 

ATTN  RDRL  CIM  G  T  LANDFRIED 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 


No.  of 

Copies  Organization 

8  US  ARMY  RSRCH  LAB 

ATTN  IMNE  ALC  HRR  MAIL  & 
RECORDS  MGMT 
ATTN  RDRL  CIM  L  TECHL  LIB 
ATTN  RDRL  CIM  P  TECHL  PUB 
ATTN  RDRL  SER  E  F  CROWN 
ATTN  RDRL  SER  U  C  FAZI 
ATTN  RDRL  SEE  M  M  WRABACK 
ATTN  RDRL  SER  A  HUNG 
ATTN  RDRL  SEEM  S  RUDIN 
ADELPHI  MD  20783-1197 

TOTAL:  16(1  ELEC,  1  CD,  14  HCS) 
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